Extended Thomas-Fermi Density Functional 
for the Unitary Fermi Gas 



00 

o 
o 

(N 

> 

o 

O 

oo 
(N 



I 

o 



> 
O 
(N 
00 

O 
00 

o 



X 



Luca Salasnich and Flavio Toigo 
Department of Physics "Galileo Galilei", CNISM and CNR-INFM, 
University of Padua, Via Marzolo 8, 35131 Padua, Italy 

We determine the energy density 5(3/5)neF and the gradient correction Xh? {SJn)^ / {9,mn) of the 
extended Thomas- Fermi (ETF) density functional, where n is number density and is Fermi 
energy, for a trapped two-components Fermi gas with infinite scattering length (unitary Fermi gas) 
on the basis of recent diffusion Monte Carlo (DMC) calculations [Phys. Rev. Lett. 99, 233201 
(2007)]. In particular we find that ^ = 0.455 and A = 0.13 give the best fit of the DMC data with 
an even number A'' of particles. We also study the odd-even splitting 'yN^^^huj of the ground-state 
energy for the unitary gas in a harmonic trap of frequency uj determining the constant 7. Finally 
we investigate the effect of the gradient term in the time-dependent ETF model by introducing 
generalized Galilei-invariant hydrodynamics equations. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10.Ay, 67.85. Lm 



I. INTRODUCTION 

The crossover from the weakly paired Bardeen-Cooper- 
SchriefFer (BCS) state to the Bose-Einstein condensate 
(BEC) of molecular dimers with ultra-cold two-hyperfine- 
componcnts Fermi vapors atoms has been investigated by 
several experimental groups with atoms 0,0; 111 and 
^Li atoms 0, IH- In the unitary limit of infinite scatter- 
ing length obtained by tuning an external background 
magnetic field near a Feshbach resonance J3, the Fermi 
superfluid exhibits universal properties 0, H M, ■ 

It has been suggested that at zero-temperature the uni- 
tary Fermi gas can be described by the density functional 
theory (DFT) 0, [H, M MmMMM- Bulgac 
and Yu [ll| introduced a superfluid DFT (SDFT) based 
on a Bogoliubov-de Gennes (BdG) approach to super- 
fluid fermions, in the same spirit as the density func- 
tional formulation for superconductors [l^. Papenbrock 
and Bhattacharyya [l3l have instead proposed a Kohn- 
Sham (KS) density functional with an effective mass to 
take into account nonlocality effects. To treat nonuni- 
form systems, other authors [H, [H, [H, [13, [III have 
added a gradient term to the leading Thomas-Fermi en- 
ergy, since such a term is surely necessary when surfaces 
are present [20I . [2l| , at least in three spatial dimensions 
p^ . An energy functional for fermions not written in 
terms of single-particle orbitals, but only in terms of 
the density and its derivatives is usually called extended 
Thomas- Fermi (ETF) functional [Ullll. It may be seen 
as an effective field theory where the gradient correc- 
tion \h'^{Vn)ll(8mn) can be interpreted as the next-to- 
leading term [la, [3i with n(r) the local number density 
and m the atomic mass. 

We wish to point out that both the energy function- 
als proposed by Bulgac and Yu |lTj and Papenbrock and 
Bhattacharyya [U are functionals of the density through 
single particle orbitals (the BdG or KS orbitals). There- 
fore they can be used in actual numerical calculations 
only when the number of fermions is small, since they re- 



quire a self consistent calculation of single-particle states 
whose number increases linearly with the number of par- 
ticles. 

On the contrary, one encounters no limitation in the 
number of particles which may be treated through ETF 
functionals, since in this case the functional depends only 
on a single function of the coordinate, i.e. the particle 
density. Of course one trades simplicity with accuracy: 
while the BdG and KS schemes are built to account for 
the main contribution to the kinetic energy, and treat 
it exactly in noninteracting systems even with a non- 
uniform density varying in space, the TF approach gives 
the exact kinetic energy only for a uniform system and 
even when extended with the addition of gradient and 
higher order derivatives of the density, the ETF func- 
tional is not able to reproduce shell effects in the density 
profile [2l,[23. 

In spite of this limitation, but in the light of the great 
simplification introduced in numerical calculations, we 
believe that it is useful to analyze the ETF approxi- 
mation and comment on its dynamical generalization, 
which amounts to introducing a quantum pressure term 
— Xh^V^^/n/ {2my/n) into the hydrodynamic equations of 
superfiuids. 

The value of the coefficient A is debated. In the papers 
of Kim and Zubarev [T^] and Manini and Salasnich (ll) 
the authors set A = 1 over the full BCS-BEC crossover. 
More recently we have suggested A = 1/4 [13, l25|] . This 
suggestion is in good agreement with a theoretical esti- 
mate based on an epsilon expansion around d = 4 — e 
spatial dimensions in the unitary regime [ist . 

In this paper we comment on the ETF for a two- 
components Fermi gas at unitarity and determine its pa- 
rameters by fitting recent Monte Carlo results [2^, [23| for 
the energy of fermions confined in a spherical harmonic 
trap of frequency lo in this regime. 

Since the interaction potential does not introduce any 
new length, the universal contribution to the energy den- 
sity in an ETF functional appropriate to a spin balanced 
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Fermi liquid at unitarity may be considered, in its sim- 
plest form, as the sum of a term proportional being 
the constant of proportionality) to the energy density of 
a uniform noninteracting system with the same density 
and of a term containing the gradient of the density with 
a coefRcient A meant to take into account phenomenolog- 
ically also higher order derivatives [1^. 

By minimizing such ETF for a fixed number of parti- 
cles, we find that the values ^ = 0.455 and A = 0.13 give 
the closest results to the Monte Carlo energies of a fully 
supcrfluid system with an even number TV of fermions 
confined in a harmonic well at unitarity, as calculated by 
Ref. [13] ■ When treating systems with an odd number of 
particles, we must correct the calculated value of the ETF 
ground state energies corresponding to these parameters 
to account for the presence of the unpaired particle. Ac- 
cording to Son [29j , for fermions confined by a harmonic 
potential, the correction depends on the number of parti- 
cle and takes the form /S.E = jN^^^ (in units of huj). We 
then find that 7 = 0.856 provides the best fit to the DMC 
data [13]. In section IV, we investigate the effect of the 
gradient term on the dynamics of the Fermi superfiuid by 
introducing generalized hydrodynamics equations and a 
Galilei-invariant nonlinear Schrodinger equation of the 
Guerra-Pusterla type [13, [Ml which is fully equivalent to 
them. 



II. EXTENDED THOMAS-FERMI 
FUNCTIONAL 

Let us consider an interacting Fermi gas trapped by a 
potential U{r). Its TF energy functional is: 



TF 



(fr n{r) e{n{r)) + U{r) 



(1) 



where e{n) is the energy per particle of a uniform Fermi 
system with density n equal to the local density n{r) of 
fermions. The total number of fermions is 



N 



d r n(r) 



(2) 



By minimizing Etf with respect to the density n(r), with 
the constraint of a fixed number of particles, one finds 



M(n(r)) + C/(r) - /I 



(3) 



where ^{n) = ^''"g^"'*'' is the bulk chemical potential of a 
uniform system of density n and p. is the chemical poten- 
tial of the non uniform system, i.e. the Lagrange multi- 
plier fixed by the normalization 

In the unitary limit, no characteristic length is set 
by the interatomic-potential since its s-wave scattering 
length ap diverges (of ^ ±00). The energy per particle 
of a uniform two-spin components Fermi gas at unitarity 
must then depend only on h, on the mass of fermions 
m, and on the only length characterizing the system, i.e. 



the average distance between particles oc n It is 

usually written as: 



s{n) 



(4) 



where ^ is a universal parameter which can be de- 
termined from ab-initio calculations. Notice that: 
(7i^/2m)(37r^)^/^n^/^ = ep, where ep is the Fermi en- 
ergy of the ideal fermionic gas. Thus ^ is simply the ratio 
between the energy per particle of the uniform interact- 
ing system at unitarity and the corresponding energy in 
a non-interacting system. Monte Carlo calculations for 
a uniform unpolarized two-spin components Fermi gas 
suggest ^ ~ 0.45; in particular, ^ = 0.42 ±0.01 according 
to g and ^ = 0.44 ± 0.02 according to The bulk 

chemical potential associated to Eq. ^ is 



2m 



(5) 



If the system is confined by a spherically-symmetric 
harmonic potential : 



C/(r) 



1 



2 2 
muj r 



its density profile n(r) obtained from Eq. ([3]) is : 

^2 \ 3/2 



n{r) = 71(0) 1 - — 



(6) 



(7) 



where n(0) = {2mflf/'^/{3TT^n^^^/'^), rp = ^2ji/{mu'^) 
and jl = hujy^{3Ny/^. Obviously, the expression ([7]) 
for the TF density profile of the unitary Fermi gas in a 
harmonic potential coincides with that of an ideal Fermi 
gas [32| , but its parameters are modified by the presence 
of ^ in the equation of state ([4]). 

As previously stressed the TF functional must be ex- 
tended to take into account other characteristic lengths 
related to the spatial variations of the density, besides 
the average particle separation. As a consequence, the 
energy per particle must contain additional terms, which 
scale as the square of the inverse of these various lengths. 
For this reason, as a simple approximation, we add to the 
energy per particle of Eq. a term 



2m V 



(8) 



which may be seen as the first term in a gradient expan- 
sion. We notice that, according to the Kirzhnits expan- 
sion of the quantum kinetic operator in powers of h [2l| , 
A must take the value A = 1/9 [H, 34] for an ideal, non- 
interacting, Fermi gas. Historically, a term of this form, 
but with A = 1, was introduced in a pioneering paper 
[20t by von Weizsacker to treat surface effects in nuclei. 

Instead than trying to determine it from first prin- 
ciples, we consider A as a phenomenological parameter 
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accounting for the increase of kinetic energy due the en- 
tire spatial variation of the density. We remark that this 
attitude is adopted in many appHcations of Density Func- 
tional Theory to atomic and molecular physics where the 
gradient term ([8]) takes into account phenomenologically, 
through A, all the possible corrections of a gradient ex- 
pansion |28| . 

The new energy functional reads [2l|, [sl, [s^, |3^ 



E 



d^r n{r) eg{n{r),Vn{r)) + U{r) 



where 



Vn) ~ e{n) + A 



(9) 



(10) 



8m "n? 

is a generalized energy per particle which includes the 
?i-dcpendcnt gradient correction. 

We remark that in our ETF the constants ^ and A 
are independent, implying that the ratios between the 
energies per particle corresponding to the two terms in a 
unitary and in a nonininteracting fermi systems may be 
different. Equal ratios imply A = f2]| . 

By minimizing the energy functional ([9]) with the 
constraint ^ one gets the partial differential equation 
obeyed by the ground state density: 

X^V^+Kn{v)) + U{v) 



^j^^nVW)- (11) 



For the study of hydrodynamics in Fermi supcrfluids 
Kim and Zubarev [l^l, and also Manini and Salasnich 
Pt . used A = 1 over the full BCS-BEC crossover. More 
recently we have suggested A = 1/4 [l3l on the basis 
of the correct relationiship between phase and superfluid 
velocity [2^. This suggestion is in agreement with the 
theoretical estimate of Rupak and Schafer [l^ , obtained 
from an epsilon expansion around d = 4— £ spatial dimen- 
sions for a Fermi gas in the unitary regime. In particular, 
by using the expansion of Rupak and Schafer [1^ we ob- 
tain the ETF functional ® with ^ = 0.475 and A = 0.25. 

Very recently two theoretical groups [1^, [l^l have stud- 
ied the two-components Fermi gas in the harmonic trap 
of Eq. ([6]) at unitarity by using Monte Carlo algo- 
rithms. Chang and Bertsch [2^ have used a Green- 
fimction Monte Carlo (GFMC) method, while Blume, 
von Stechcr and Greene [l^l have applied a fixed-node 
diffusion Monte Carlo (FN-DMC) approach. They have 
obtained the ground-state energy for increasing values of 
the total number N of fermions. 

Notice that by inserting the Thomas-Fermi profile ([7]) 
into Eq. ([9]) one gets the ground-state energy of N 
fermions in a harmonic potential with frequency lo in the 
form: 

_E _ ^ {mf/^ 9A(37V)2/3~ 

We will refer to this expression as the "beyond-TF" en- 
ergy. The first term is the TF contribution to the ground- 
state energy, while the second term is the leading correc- 
tion. Chang and Bertsch [2^ and Blume, von Stechcr 



(12) 




FIG. 1: (Color online). Unitary Fermi gas under harmonic 
confinement of frequency uj. Upper panel: density profiles 
n{r) for iV = 10 and iV = 30 fermions obtained with ETF 
(solid lines) and TF (dashed lines). Lower panel: ground- 
state energy E vs N with ETF (solid line), beyond-TF for- 
mula (dot-dashed line) and TF (dashed line). In all calcu- 
lations: universal parameter ^ = 0.44 and gradient coeffi- 
cient A — 1/4. Energy in units of hu! and lengths in units of 
an = y/h/{muj). 



and Greene [27| have determined ^ by fitting their MC 
data with Eq. p2)) under the hypothesis that the relation 



A 



(13) 



appropriate to a noninteracting system holds also in 
the unitary regime. They find respectively the values 
^ = 0.50 III] and £, = 0.465 [13, which are compatible 
with previous determinations of the parameter ^ based on 
Monte Carlo calculations for the uniform svstemjgl. [l0|. 
The corresponding value A = ^/9 ~ 0.05 for the coeffi- 
cient of the gradient is instead much smaller than previ- 
ous suggestions [H, El, [13, [H, HI . 

It is important to stress that the formula in Eq. (fT^ 
does not give the minimum of the ETF functional 
since it corresponds to the density profile of Eq. ^ and 
not to the true ground-state density, solution of Eq. (|lip . 

In the upper panel of Fig. 1 we plot ETF density pro- 
files (solid lines) and compare them with the TF ones 
(dashed lines). The ETF density profiles have been de- 
termined by solving Eq. pT|) with ^ = 0.44 and A = 1/4 
with a finite-difference numerical code [36j As expected, 
there are visible differences, in particular near the sur- 
face. In the lower panel of Fig. 1 we plot the ground- 
state energy E for increasing values of the number N of 
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FIG. 2: (Color online). Ground-state energy E for the uni- 
tary Fermi gas of A'^ atoms under harmonic confinement of 
frequency u). Symbols: DMC data with even A'' [SJ; solid 
line: ETF results with best fit = 0.455 and A = 0.13); 
dot-dashed line: ETF results obtained from e-expansion jlSj l 
= 0.475 and A — 0.25). Energy in units of hco. 



fermions. Here the differences between TF (dashed line), 
beyond-TF (dot-dashed line) and ETF (solid hne) are 
quite large. The figure clearly shows that the beyond-TF 
formula (fT2|) is not very accurate. 

We stress here once again that the values of the pa- 
rameters ^ and A in the ETF functional should be uni- 
versal i.e. independent on the confining potential U (r) 
[H, nil . Moreover we consider A as taking into account 
phenomenologically all possible corrections of a gradi- 
ent expansion in the unitary regime and treat ^ and A 
as independent parameters. To determine them, instead 
than using the inaccurate Eq. (fT2|) , we look for the values 
of the two parameters which lead to the best fit of the 
FN-DMC ground-state energies [13] for even N. After a 
systematic analysis we find ^ = 0.455 and A = 0.13 as the 
best parameters in the unitary regime. It is important to 
observe that the value ^ = 0.455 of our best fit coincides 
with that obtained by Perali, Fieri and Strinati [ssl ] by 
using the extended BCS theory with beyond-mean-field 
pairing fluctuations. 

In Fig. 2 we plot the ground-state energy E of the 
Fermi gas under harmonic confinement, comparing dif- 
ferent results: the FN-DMC data for an even number N 
of atoms of [l^l (diamonds with error bar), the best ETF 
results, with A = 0.13 and ^ = 0.455 (solid line), and the 
ETF results obtained by using the values £ = 0.475 and 
A = 0.25 coming from the e-expansion [18[ (dot-dashed 
line) . 

We remark that fixing £ to the value £ — 0.44 and 
looking for the best fitting A we have found A = 0.18. 
In this case the curve of the energy will be practically 
superimposed to the solid one of Fig. 2. 

For the sake of completeness in Table 1 wc report the 



fixed node DMC energies of Blume, von Stecher and C. 
H. Greene our optimized ETF results with £ = 0.575 
and A = 0.13, and also the SDFT calculations of Bulgac 
[U. Remarkably the ETF energies are slightly closer to 
the DMC values than the SDFT ones, reported in Ref. 
[stI ]. The optimized ETF energies for even number N 
of particles are obtained from our density functional ([9|. 
while the energies with odd number N of particles are 
calculated taking into account the odd-even splitting, as 
discussed in the next section. 



TV 


Edmc 


Eetf 


EsDFT 


2 


2.002 


2.17 


2.33 


3 


4.281 


4.65 


4.62 


4 


5.051 


5.19 


5.52 


5 


7.61 


7.98 


7.98 


6 


8.64 


8.71 


9.07 


7 


11.36 


11.73 


11.83 


8 


12.58 


12.61 


12.94 


9 


15.69 


15.81 


16.06 


10 


16.80 


16.83 


17.15 


11 


20.11 


20.19 


20.36 


12 


21.28 


21.32 


21.63 


13 


24.79 


24.82 


24.96 


14 


25.92 


26.04 


26.32 


15 


29.59 


29.67 


29.78 


16 


30.88 


30.99 


31.21 


17 


34.64 


34.73 


34.81 


18 


35.96 


36.13 


36.27 


19 


39.83 


39.99 


40.02 


20 


41.30 


41.46 


41.51 


21 


45.47 


45.41 


45.42 


22 


46.89 


46.96 


46.92 


23 


51.01 


51.01 




24 


52.62 


52.63 




25 


56.85 


56.76 




26 


58.55 


58.45 




27 


63.24 


62.66 




28 


64.39 


64.41 




29 


69.13 


68.70 




30 


70.93 


70.51 





Table 1. Ground-state energy E of the unitary Fermi gas 
of A*' even atoms under harmonic confinement, in units of 
the harmonic energy Tiuj. Comparison among different cal- 
culations: fixed node diffusion Monte Carlo [53] {Edmc), 
our optimized extended Thomas-Fermi with ^ = 0.455 and 
A = 0.13 (Eetf), and the superfluid density functional the- 
ory [13 {Esdft). 
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FIG. 3: (Color online). Ground-state energy E for the uni- 
tary Fermi gas of A'^ atoms under harmonic confinement of 
frequency u. Dimonds: DMC data with both even and odd 
N [23; solid line: optimized ETF results = 0.455, A = 0.13, 
7 = 0.856). Energy in units of hui. 



III. ODD-EVEN SPLITTING 

Up to now wc have analyzed the unitary gas with an 
even number N of particles. The Monte Carlo calcula- 
tions [2^, [131 show a clear odd-even effect, reminiscent of 
the behavior of the nuclear binding energy. In particular, 
denoting by E]\f the ground state energy of N particles 
in a isotropic harmonic trap, for odd N one finds 



En — -^{En-i + En+i) + Atv 



(14) 



where the splitting Ajv is always positive. This effect is 
related to pairing: given the superfluid cloud of even par- 
ticles, the extra particle is localized where the energy gap 
is smallest, which is near the edge of the cloud [2^, [37| . 
On this basis, recently Son [29| has suggested that, for 
fermions at unitarity, confined by a harmonic potential 
with frequency w, the odd-even splitting grows as 



AEn = 7 N^^^ TiLo , 



(15) 



where 7 is an unknown dimensionless constant. After a 
systematic investigation we find that 7 — 0.856 gives the 
odd-even splitting which best fit entire FN-DMC data, 
with both even and odd particles [IJ. In Fig. 3 we report 
the FN-DMC data (diamonds) and the optimized ETF 
results (solid line). The figure, which displays the zig-zag 
behavior of the energy as a function of iV, shows that 
the optimized ETF functional plus the odd-even correc- 
tion (Hil) (C = 0.455, A = 0.13, 7 = 0.856) is extremely 
good in reproducing all FN-DMC data. 



IV. GENERALIZED HYDRODYNAMICS 

Let us now analyze the effect of the gradient term 
on the dynamics of the unitary Fermi gas. At zero 
temperature the low-energy collective dynamics of this 
fermionic gas can be described by the equations of gen- 
eralized hydrodynamics [3^, [H, [l^l , where the Hamilto- 
nian of the classical hydrodynamics [40[ is modified by 
including gradient corrections. In our case the general- 
ized Hamiltonian reads: 



H 



(TV n 



-mv^ + eg{n, Vn) + U (r) 



(16) 



where the local density n{r,t) and the local velocity 
v(r, t) are the hydrodynamics variables [H, [10]. By writ- 
ing the Poisson brackets [13 of the hydrodynamics vari- 
ables with the Hamiltonian p6)) . one gets the generalized 
hydrodynamics equations: 



dn 

In 



9 

m|-+v.V 



V 



V • (nv; 
tig{n,Wn) + U{r) 



0, 

= 



(17) 
(18) 



where 



d{neg{n,Vn)) 
fig{n, Vn) = — = ^i(n) - A 



2m J\ 



(19) 

Eq. (fT7|) is the continuity equation, while Eq. is the 
conservation of linear momentum. These equations are 
valid for the inviscid unitary Fermi gas at zero tempera- 
ture. If the unitary Fermi gas is superfiuid then it is not 
only inviscid but it is also irrotational, i.e. 



V X V = . 
By using this condition and the identity 

(v • V)v = V i^,""^ - V X (V X v) 

Eq. can be simplified into: 
■f 



(20) 



(21) 



9v „ 

at 



-mw^ + /ig (n, Vn) + U (r) 



. (22) 



The low-energy collective dynamics of the superfluid 
Fermi gas in the BCS-BEC crossover is usually described 
by the equations of classical superfluid hydrodynamics, 
which are the time-dependent version of the local den- 
sity approximation with the Thomas-Fcrmi energy func- 
tional [3]. Equations (fTTjl and (|22|) arc a simple general- 
ization of classical superfluid hydrodynamics which takes 
into account surface effects. The gradient term, i.e. the 
quantum correction, is necessary in a superfluid to avoid 
unphysical phenomena like the formation of wave front 
sing ularities in the dynamics of dispersive shock waves 

s 
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Combining Eqs. ((T7)) and (|22p one finds the dispersion 
relation of low-energy collective modes of the uniform 
unitary Fermi gas in the form 



Finally, we stress that Eq. ((27|) can be seen as 
the Euler-Lagrange equation of the following Galilei- 
invariant Lagrangian density 




1 + 

g V2toco 



(23) 



where is the collective frequency, q is the wave number 
and Co is the speed of sound in a uniform, noninteracting 
Fermi gas. 

For an irrotational fluid it is possible to write down 
[iol . a Lagrangian by using as dynamical variables 
the scalar potential 0{r,t) of the velocity v(r,<) and the 
local density n(r,t). For a fermionic superfluid one has: 



-£,(|»P,V(|*p))|*|» + — (V|«|)^ (28) 

that is equivalent to the generalized Lagrangian of Eq. 
(El). 



v(r,t) = -V.(r,t) 



(24) 



where 0{r,t) is the phase of the condensate wavefunc- 
tion of Cooper pairs 0. In our case, the familiar La- 
grangian of the Fermi superfluid [H, \^ must be modi- 
fied by including the gradient correction. The generalized 
Lagrangian density then reads: 

£ = _„g0+^(V0)2 + c/(r)+e,(n,Vn)) . (25) 

The Euler-Lagrange equations of this Lagrangian with 
respect of n and 9 give the generalized hydrodynamics 
equations of superfiuids (fT7|) and ((22)) . 

We observe that the generalized hydrodynamics equa- 
tions (fT7|) and can be formally written in terms of 
a nonlinear Schrodinger equation of the Guerra-Pustcrla 
type [soj . which is Galilei-invariant [3l|. In fact, by in- 
troducing the complex wave function 



*(r,t) = V^^Mje'^^"-'*) , 



(26) 



which is the zero-temperature Ginzburg-Landau order 
parameter normalized to the total number N of super- 
fluid atoms HI] , and taking into account the correct 
phase- velocity relationship given by Eq. ([M]) . the equa- 
tion 



= 



4m 



22//(|1'n + (l-4A) 



4m 1*1 



* , (27) 



is strictly equivalent to Eqs. (fT7|) and (|22|) . Notice that 
in the stationary case where 5'(r, i) = ^ n{r) e^'-'^ii^/f^ ^ 
Eq. l[27|) becomes exactly Eq. (fTT|) . Remarkably only if 
A = 1/4 the equation acquires the familiar structure of a 
nonlinear Schrodinger such as the Gross-Pitaevskii equa- 
tion which describes the two-spin component Fermi sys- 
tem, but in the extreme BEG regime [4j| . From the lin- 
earization of Eq. (|27p one finds for the uniform Fermi gas 
the Bogoliubov excitations given precisely by Eq. 



V. CONCLUSIONS 

We have obtained the value of the coefficient A of the 
gradient correction Xh^ {V n)^ / (Smn) for the extended 
Thomas-Fermi density functional in the unitary regime. 
By fitting diffusion Monte Garlo data with an even num- 
ber TV of particles we have found A = 0.13. In addition, 
we have determined the coefficient g of the energy den- 
sity g(3/5)nei?, finding g = 0.455. Fixing g to the value 
g = 0.44 proposed in [l^l and looking for the best fitting A 
we have found instead A = 0.18. We stress that in our en- 
ergy functional the gradient term takes into account phe- 
nomenologically all corrections of a gradient expansion. 
Our functional one can easily get the ground state prop- 
erties (energy and density) for large as for small numbers 
of fcrmions; its main limitation is that it cannot account 
for the shell effects in the density profile. Moreover, we 
have shown that it is possible to take into account the 
odd-even splitting of the ground-state energy of the uni- 
tary gas in a harmonic trap of frequency uj by considering 
a correction proportional to N^/^huj as suggested by Son 
[29| . where the constant of proportionality is found to 
be 7 = 0.856. Finally, we have analyzed the effect of 
the gradient term in the dynamics of the unitary Fermi 
gas by introducing generalized hydrodynamics equations, 
which can be written for superfiuid motion in the form 
of a Galilei-invariant nonlinear Schrodinger equation. As 
a final remark, we remind that the values of g and A 
we have found are independent on the external potential 
and therefore our generalized energy functional can be 
used to investigate the unitary superfluid Fermi gas in 
various trapping configurations. Obviously in the limit 
of large numbers of particles the gradient term and the 
exact value of A are less important since the dominant 
term becomes the usual Thomas-Fermi one. More exten- 
sive Monte-Carlo calculations with a larger number of 
particles will be certainly useful to have a more accurate 
determination of the value of g. 
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